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The linear sigma model with quarks at very low temperatures provides an effective description 
for the thermodynamics of the strong interaction in cold and dense matter, being especially useful 
at densities found in compact stars and protoneutron star matter. Using the MS one-loop effective 
potential, we compute quantities that are relevant in the process of nucleation of droplets of quark 
matter in this scenario. In particular, we show that the model predicts a surface tension of E ~ 5- 
15 MeV/fm'^, rendering nucleation of quark matter possible during the early post-bounce stage of 
core collapse supernovae. Including temperature effects and vacuum logarithmic corrections, we 
find a clear competition between these features in characterizing the dynamics of the chiral phase 
conversion, so that if the temperature is low enough the consistent inclusion of vacuum corrections 
could help preventing the nucleation of quark matter during the collapse process. We also discuss 
the first interaction corrections that come about at two-loop order. 



I. INTRODUCTION 



The thermodynamics of strong interactions for cold 
matter under extremely high densities is of utmost im- 
portance for the understanding of the structure of com- 
pact stars Since quantum chromodynamics (QCD) 
is asymptotically free, it is believed that for high enough 
densities quarks will be in a deconfined state, the quark- 
gluon plasma |2j. Moreover, due to the approximately 
chiral nature of the QCD action, one also expects quarks 
to be essentially massless above a sufHciently high value 
of the chemical potential^. Depending on the location of 
the critical density, one might find several sorts of con- 
densates and even deconfined quark matter in the core of 
neutron stars ^Sj. Furthermore, the order and strength 
of the chiral transition are crucial features in establishing 
the existence of a new class of compact stars @, 0] ■ 

Recently, it was shown that deconfinement can happen 
early, during the early post-bounce accretion stage of a 
core collapse supernova event, which could result not only 
in a delayed explosion but also in a neutrino burst that 
could provide a signal of the presence of quark matter in 
compact stars [s^]. However, as was discussed in detail 
in Ref. 9] (see also Ref. [l^) those possibilities de- 
pend on the actual dynamics of phase conversion, more 
specifically on the time scales that emerge. In a first- 
order phase transition, as is expected to be the case in 
QCD at very low temperatures, the process is guided by 
bubble nucleation (usually slow) or spinodal decomposi- 
tion ("explosive" due to the vanishing barrier), depend- 
ing on how fast the system reaches the spinodal insta- 
bility as compared to the nucleation rate. Nucleation in 
relatively high-density, cold strongly interacting matter. 
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^ Nevertheless, these features and other thermodynamic properties 
of dense media were also shown to be significantly affected by 
nonzero quark masses . 



with chemical potential of the order of the temperature, 
can also play an important role in the scenario proposed 
in Ref. [11], where a second (little) infiation at the time 
of the primordial quark-hadron transition could account 
for the dilution of an initially high ratio of baryon to 
photon numbers. A key ingredient in both scenarios is, 
of course, the surface tension, since it represents the price 
in energy one has to pay for the mere existence of a given 
bubble (or droplet). 

However, the surface tension and the whole process of 
phase conversion via nucleation in a cold and dense envi- 
ronment are not known for strong interactions. In fact, 
the mapping of this sector of the (equilibrium) phase di- 
agram is still in its infancy (l^ , let alone dynamical pro- 
cesses of phase conversion. This region of the phase dia- 
gram is not very amenable to first-principle approaches. 
On one hand, the so-called sign problem at finite chemi- 
cal potential brings about major technical difficulties for 
performing Monte Carlo lattice simulations, so that one 
still does not have the high-quality guidance from lat- 
tice QCD that is currently available at finite tempera- 
ture and zero density [l3]- On the other hand, although 
the perturbative series for QCD at zero temperature and 
finite density seems to be much better behaved than its 
counterpart at finite temperature ,§,,14.- 16], the values of 
chemical potential that are phenomenologically interest- 
ing in the interior of compact stars are already too low, of 
the order of 400 — 500 MeV 1]. At this scale, one is cer- 
tainly pushing perturbative QCD much beyond its limits 
of applicability. Besides, this approach is incapable of in- 
corporating the chiral and the deconfinement transitions, 
as well as vacuum properties, unless it is merged with a 
complementary low-density treatment. Therefore, an es- 
timate of the surface tension and the description of the 
nucleation process within an effective model for strong 
interactions in a cold and dense environment seems wel- 
come. 

In this paper we consider the linear sigma model cou- 
pled with constituent quarks (LSMq) at zero or low tem- 
perature and finite quark chemical potential as a model 
for the thermodynamics of strong interactions in cold and 
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dense matter. We compute the effective action for the 
sigma condensate integrating over the quark fields and 
keeping quadratic fluctuations of the sigma field around 
the condensate. Our final thermodynamic potential in- 
corporates all corrections from the medium and vacuum 
fluctuations in the MS scheme, including logarithmic 
and scale-dependent contributions from quark and sigma 
bubble-diagrams. Renormalization is implemented in the 
standard fashion in the MS scheme. Knowing the full 
thermodynamic potential, we study the behavior of the 
chiral condensate with chemical potential, the location 
and nature of the chiral transition, and the process of 
homogeneous nucleation. 

The paper is organized as follows. Section II presents 
the LSMq and some of its features, and the analytic com- 
putation of the one- loop effective potential at flnite chem- 
ical potential, including quantum logarithmic corrections 
in the vaccum and thermal effects. In Section III we 
show our results for the surface tension and other rele- 
vant quantities for the nucleation process, analyzing the 
competition between vacuum and thermal modifications 
and discussing the consequences for the supernovae ex- 
plosion scenario. Section IV contains our conclusions and 
outlook. 



II. EFFECTIVE THEORY 
A. General framework 

The LSMq, also known as the quark-meson model, is 
very suitable for the study of the chiral transition [17|. 
As argued in Ref. [18|, QCD with two fiavors of mass- 
less quarks belongs to the same universality class as the 
0(4) LSM, exhibiting the same qualitative behavior at 
criticality. Besides, the LSM is renormalizable [l^ and 
reproduces correctly the phenomenology of QCD at low 
energies, such as the spontaneous (and small explicit) 
breaking of chiral symmetry, meson masses, etc. In spite 
of the fact that an effective theory does not require renor- 
malizability to be well posed, this attribute is highly de- 
sirable if one is interested in investigating the behavior of 
physical quantities as the energy scale is modified, which 
can be accomplished via renormalization group methods. 

Since its proposal [l^, the LSM has been investigated 
in different contexts, from the low-energy nuclear theory 
of nucleon-meson interactions to ultra-relativistic high- 
energy heavy-ion collisions, and also in different varieties, 
e.g. including explicitly constituent quark degrees of free- 
dom or not. The thermodynamical aspects of the model 
were first considered in the very early days of finite- 
temperature field theory [2l|, and systematic approxi- 
mations for the study of the chiral transition started to 
be implemented soon afterwards (2^ . producing an ex- 
tensive literature. 

When used to mimic and study the chiral transition of 
QCD, the emphasis was generally on thermal effects [22I - 
[29j| rather than considering a cold and dense scenario, al- 



though chiral symmetry restoration at high densities was 
predicted quite early |30j. This choice was, of course, well 
justified by the stimulating experimental results coming 
from high-energy heavy-ion collisions [3ll|, and by the 
possibility to compare to numerical output from lattice 
QCD [1^. Usually the gas of quarks is treated as a ther- 
mal bath in which the long- wavelength modes of the chi- 
ral field evolve, the latter playing the role of an order pa- 
rameter in a Landau-Ginzburg approach. The standard 
procedure is then integrating over the fermionic degrees 
of freedom, using a classical approximation for the chi- 
ral field, to obtain a formal expression for the thermody- 
namic potential from which one can compute all the phys- 
ical quantities of interest. In this case, the sigma field 
is approximated by the condensate, and the functional 
integral over sigma fluctuations is not performed. The 
fermionic determinant is usually calculated to one-loop 
order assuming a homogeneous and static background 
field m. 

In a theory with spontaneous symmetry breaking, the 
presence of a condensate will modify the masses. In par- 
ticular, in the case of the LSMq the masses of quarks and 
mesons will incorporate corrections that are functions 
of the chiral condensate, which is a medium-dependent 
quantity. Therefore, contributions from vacuum dia- 
grams can not be subtracted away as trivial zero-point 
energies since they contain medium-dependent pieces 
via effective masses. So, although the presence of the 
medium brings no new ultraviolet divergence, in principle 
one has to incorporate carefully finite vacuum contribu- 
tions that survive the renormalization procedure. These 
contributions have been very often discarded in studies 
of the LSMq, but were shown to play an important role 
by the authors of Ref. [s^ who incorporate scale effects 
phenomenologically. Vacuum contributions were consid- 
ered at finite density in the perturbative massive Yukawa 
model with analytic exact results up to two loops in Ref. 
34, ,35,] and, more specifically, in optimized perturbation 
theory at finite temperature and chemical potential in 
Ref. [36|, also comparing to mean-field theory. More re- 
cently, this issue was discussed in a comparison with the 
Nambu-Jona-Lasinio model [s^l, in the Polyakov-LSMq 
model in the presence of an external magnetic field [ssf . 
and in the quark-meson model, with special attention to 
the chiral limit fsot- 

To describe the chiral phase structure of strong inter- 
actions at finite density and vanishing temperature we 
adopt the LSMq, defined by the following lagrangian 

£ = ipf [i^'^d^ - ruf - g{cr + i-f^f ■ Tf)]->pf 

+ ^{d^ad^a + d^7r-d''Tf)-U{a,7T) , (1) 

where 

U{a,Tf) = ^{a^ +Tr^ -v^f -ha (2) 

is the self-interaction potential for the mesons, exhibiting 
both spontaneous and explicit breaking of chiral symme- 
try. The Nf = 2 massive fermion fields ipf represent 
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the up and down constituent-quark fields ■0/ = {u,d). 
For simplicity, we attribute the same mass, mf — rUg, 
to both quarks. The scalar field a plays the role of an 
approximate order parameter for the chiral transition, 
being an exact order parameter for massless quarks and 
pions. The latter are represented by the pseudoscalar 
field TT = (tti, 7r2, TTa). It is customary to group together 
these meson fields into a 0(4) chiral field cj) = (cr, tt). 
For simplicity, we discard the pion dynamics from the 
outset, knowing that they do not affect appreciably the 
phase conversion process [23| , and focus our discussion in 
the quark-sigma sector of the theory. Nevertheless, pion 
vacuum properties will be needed to fix the parameters 
of the lagrangian later. As will be discussed below, the 
parameters of the lagrangian are chosen such that the 
effective model reproduces correctly the phenomenology 
of QCD at low energies and in the vacuum, such as the 
spontaneous (and small explicit) breaking of chiral sym- 
metry and experimentally measured meson masses. 

Since chiral symmetry is spontaneously broken in the 
vacuum, it is convenient to expand the a field around the 
condensate, writing a{x, r) = {a)+^{x, t), where r is the 
imaginary time in the Matsubara finite-temperature for- 
malism. Although we assume the system to be at tem- 
perature T = 0, it is convenient to compute the loop 
expansion at finite temperature, taking the limit T — > 
at the end. For a dense system, the chiral condensate, 
(cr), will be a function of the quark chemical potential, /i. 
Given the shift above, the fluctuation field ^ is such that 
{£,) = and ^(p = 0,a; = 0) = 0. From the phenomenol- 
ogy, one expects that (cr)(/i oo) (being equal in 
the case of massless quarks). 

Keeping terms up to O(C^), one obtains the following 
effective lagrangian 

+ Id^^d^^-^M^'^Uiia)), (3) 
where we have defined the /i-dependent effective masses 



standard methods of finite-temperature field theory 



B. Cold and dense 1-loop effective potential and 
vacuum corrections 

The effective potential can be computed exactly and 
in closed form following the steps detailed in Ref. |3j], 
a procedure that can be even generalized to optimized 
perturbation theory (36| . Keeping only contributions to 
one loop order, the effective potential is given by 



(5) 



where r2™"(tT) is the fully-renormalized thermodynamic 
potential for the fluctuation effective theory. The latter 
corresponds to a Yukawa theory for massive fermions and 
a massive scalar, with masses given by ([4]) as functions of 
a. The thermodynamics of this theory was fully solved 
analytically to two loops in the cold and dense regime 
in Ref. (SJ], where all details can be found. In the MS 
scheme, fl^™{a) can be written as a sum of a medium 
contribution 
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-NfN^ 
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and a vacuum contribution 
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- Mg log ^ ^-jp- ]t Pf is the Fermi momen- 



where u = fJ^Pf- 

tum and A is the MS scale. The latter can also be fixed 
by vacuum properties, as will be discussed in subsection 
EDI 



Mg = mg+g{a) , = 3{a)^ - Xi 



(4) 



Linear terms in ^ can be dropped in the action because 
of the condition ^{p = 0,a; = 0) = 0. 

In the effective lagrangian, the medium-dependent 
masses will entangle vacuum and medium contributions 
in the loop expansion, rendering the renormalization pro- 
cess more subtle. This is a consequence of the pres- 
ence of a nonzero /i-dependent condensate in the bro- 
ken phase. Of course, the ultraviolet divergences are the 
ones of the theory in the vacuum, and renormalization is 
implemented in the standard fashion by adding medium- 
independent counterterms to the original lagrangian ([T|) . 
In our case, though, it is more convenient to intro- 
duce counterterms in the effective theory, so that there 
are contributions from pure-vacuum, pure-medium and 
vacuum-medium pieces in the cancellation of ultraviolet 
divergences. Renormalization is then implemented by us- 



C. Thermal effects 

We can also incorporate thermal effects in the calcu- 
lation of the effective potential within the LSMq. The 
inclusion of the temperature dependence allows for test- 
ing the validity of the cold (T = 0) approximation at 
low temperatures and to investigate if the thermal cor- 
rections to nucleation parameters can play an important 
role in the phase conversion process. 

At one loop, the well-known temperature- and density- 
dependent correction to the thermodynamic potential is 
that of an ideal gas of massive sigma particles and con- 
stituent quarks (w^ = -f and = p'^ + M^): 



a 
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ed,Th 



T 



(2^ 



log 
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-2TNfN, 



dh 



-(E,-p,)/T 



+ log[l + e-(^'+'')/^]} . (8) 

D. Parameter fixing 

As stated above, the LSMq is adopted as an effective 
model for QCD at low energies, so that the parameters 
5, A, TOq, w, h and A are fixed in order to reproduce QCD 
properties either measured in the vacuum or calculated 
numerically via lattice QCD. Therefore, the conditions 
for fixing the parameters are imposed on the vacuum ef- 
fective potential. Since we aim at comparing results from 
cases with different vacuum effective potentials (namely 
U and U -|-fiiac), and even with different parameter sets^, 
it is useful to make explicit the parameter fixing proce- 
dure, consistently. 

The conditions for fixing the model parameters are the 
following: 

• The chiral condensate in the vacuum is the pion 
decay constant, = 93 MeV, or, in terms of the 
minimum of the vacuum effective potential. 



d{a) 



= 0; 



(9) 



The value of the dressed mass of the a field is given 
by the experimental value of the mass of the a me- 



a2T/vac 



d{a) 



{Ml'^y « (600 MeV)^ 



(13) 



• The quark condensate is fixed by the lattice result 
(including only quarks up and down, we have [40j: 
{i^f^f)ya.c = -2 (225 MeV)3), so that 



QT/vac 



dnia 



(14) 



The conditions pil)) - (|12l) above fix directly the param- 
eters /i, ruq and g, independently of the inclusion of quan- 
tum corrections to the vacuum thermodynamic poten- 
tial. On the other hand, the Eqs. dS]), (US]) and are 
coupled equations to determine the parameters v and A 
(and A, if quantum corrections to the vacuum are con- 
sidered) and depend on the explicit form of the vacuum 
effective potential. In the case without quantum cor- 
rections, i.e. with the vacuum effective potential being 
purely the classical potential (V^™'^ — U),we find A^ ~ 20 
and v'^ ss 7696.8 MeV^. With the addition of the 1-loop 
; to the vacuum effective potential, the 
dU, (HI and dUl) yields A^ « 16.65, 



(1) 



vacuum term O 
solution of Eqs. 

« 3296.89 MeV^ and A « 16.48 MeV. 



• The partial conservation of the axial current yields 

h = U ml = (93 MeV) (138 MeV)^ , (10) 
where is the pion mass; 

• The current quark masses calculated in lattice QCD 
(cf. e.g. Ref. [40]) provide: 

nig = 4.1 MeV , (11) 

which we neglect^, setting niq = 0, for the sake of 
comparison with previous model calculations; 

• Using the constituent quark mass in the vacuum as 
1/3 of the nucleon mass (niN — 938 MeV), we can 
fix the Yukawa coupling, g: 

Af,((a)(p/=0)) = M-^ 

3=7- ( ^ rTiAT - ) = 3.32 ; (12) 



^ The case with vacuum logarithmic terms has an extra parameter: 
the MS scale A. 

^ The current quark mass also does not bring any extra qualitative 
feature to the model nor changes significantly the quantitative 
results concerning the chiral phase transition i35|] . 



E. Influence of interactions: higher-loop 
corrections 

So far, we have only included the effects of interactions 
indirectly in the construction of the effective model itself, 
through dressed masses and the presence of a nonzero 
quark condensate. However, the interaction of the sigma 
meson with the medium constituent quarks could in prin- 
ciple alter the predicted dynamics for the chiral transi- 
tion. The incorporation of such corrections in the calcu- 
lation of the effective potential can be implemented via 
the perturbative technique order by order. 

The first interaction correction in the present case ap- 
pears at the two-loop order. The contribution to the 
effective potential corresponds then to that of the ther- 
modynamic potential of an interacting Yukawa theory 
with massive scalars and massive fermions, which was ob- 
tained and analyzed in detail in Ref. [11], including the 
vacuum contribution and nonperturbative effects within 
the optimized perturbation theory framework. There- 
fore, one has in principle all the machinery to investigate 
the infiuence of interactions on the phase diagram of the 
LSMq and the associated process of homogeneous nucle- 
ation. It should be noted that the full case, including 
quantum corrections in the vacuum effective potential, 
requires the (nontrivial) solution of the coupled equa- 
tions ®, (HSl) and ((HI) with the 2- loop resuh plugged 
in. Being a more technical analysis, it is out of the scope 
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of this work, concentrated on quantifying approximately 
the nucleation predicted within an effective model for 
strong interactions at low energies. We postpone it for a 
future publication. 



III. SURFACE TENSION AND NUCLEATION 



where the minima for the symmetric and broken phases, 
as well as the barrier between them, are located. 

A quartic potential such as Eq. ([T5|) can always be 
rewritten in the form 

V[v)^a{^^ -a^f , (16) 

with the coefficients above defined as follows: 



Let us now consider the formation of droplets of quark 
matter at high density that happens via homogeneous nu- 
cleation [4l|. Dynamically this process will occur either 
via thermal activation of droplets (thermal nucleation) 
or quantum nucleation. The physical setting we have 
in mind is the one that gives the best chances for ther- 
mal nucleation of quark droplets in cold hadronic matter 
found in "hot" protoneutron stars. As discussed in Ref. 
0, and previously in Refs. d^H^, that corresponds to 
temperatures of the order of 10 — 20 MeV [43|. At these 
temperatures, and in the presence of a barrier in the effec- 
tive potential, thermal nucleation dominates over quan- 
tum nucleation. As soon as the barrier disappears, the 
spinodal instability is reached and the mechanism that 
takes over is the explosive spinodal decomposition. The 
range of temperatures under consideration is, then, high 
enough to allow for thermal nucleation (quantum nucle- 
ation being comparatively negligible 0, E3, El]) and low 
enough to justify the use of the zero-temperature effec- 
tive potential computed previously. Temperatures of a 
couple of tens of MeV will not modify appreciably the 
equation of state, bringing corrections 0(r^//^^) ~ 1% 
for the typical values of chemical potential for the system 
under consideration. Nevertheless, as will be shown later, 
thermal corrections can be important for the process of 
nucleation. 



A. Method for extracting nucleation parameters 
from the effective potential 

Since our framework is an effective model, we can only 
aim for reasonable estimates and functional behavior, not 
numerical precision. Therefore, it is convenient to work 
with approximate analytic relations by fitting the rele- 
vant region of the effective potential by a quartic poly- 
nomial and working in the thin- wall limit approximation 
for bubble nucleation. Following Refs. [23,140^, we can 
express the effective potential over the range between the 
critical chemical potential, /ic, and the spinodal, /^sp, in 
the familiar Landau-Ginzburg form 



(15) 



Although this approximation is obviously incapable of re- 
producing all three minima of 14ff, this polynomial form 
is found to provide a good quantitative description of 
this function in the region of interest for nucleation, i.e. 
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(17) 
(18) 

(19) 
(20) 



The new potential V{(p) reproduces the original Ves{4') 
up to a shift in the zero of energy. We are interested in 
the effective potential only between /x^ and ^sp- At /icj 
we will have two distinct minima of equal depth. This 
clearly corresponds to the choice j = in Eq. (|16p so 
that V has minima ai ip — ±a and a maximum at = 0. 
The minimum at ip — ~a and the maximum move closer 
together as the chemical potential is shifted and merge 
at fisp- Thus, the spinodal requires jjac? = — 8/3-\/3 
in Eq. (jl6p . The parameter j /aa^ falls roughly linearly 
from 0, at /i — fic, to — 8/3\/3 at the spinodal. 

The explicit form of the critical bubble in the thin- wall 
limit is then given by [l^l 



1 



W2a 



1 — tanh 



Rc 



(21) 



where is the new false vacuum, Rc is the radius of 
the critical bubble, and <^ — 2/m, with = V"{iff), 
is a measure of the wall thickness. The thin-wall limit 
corresponds to £,/Rc ^ 1 [ij, which can be rewritten as 
{3\j\/8aa^) < 1. Nevertheless, it was shown in j27l. lisj 
for the case of zero density and finite temperature that 
the thin-wall limit becomes very imprecise as one ap- 
proaches the spinodal (This is actually a very general 
feature of this description [4l[).In this vein, also as re- 
marked above, the analysis presented below is to be re- 
garded as semi-quantitative, and we aim for estimates. 

In terms of the parameters a, a, and j defined above, 
we find [13,141 



iftj w ±a- 



1 



1/2 



(22) 
(23) 



a(3(/?j — 

in the thin-wall limit. The surface tension, S, is given by 



E = 



dr 



dipt 
dr 



(24) 
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and the critical radius is obtained from Rc = {2T,/AV), 
where AV = V{(l)f) - V{(j)t) ~ 2a\j\. Finally, the free 
energy of a critical bubble is given by Ft = (47rS/3)i?^, 
and from knowledge of Fi, one can evaluate the nucleation 
rate F ^ e^^''^'^ . In calculating thin- wall properties, we 
shall use the approximate forms for 0t, E, and Ay 
for all values of the potential parameters. 

B. Results for nucleation in the cold and dense 
LSM 

In what follows, we consider the LSM with quarks in 
the absence of vacuum corrections. The effective poten- 
tial up to 1-loop order is then given by Eq. ([5]) with 
fi^™ = ^mod- 'This corresponds to the standard case, 
adopted frequently in the literature (cf. e.g. Ref. [26j). 

Using the method described in the previous section, 
we characterize quantitatively the nucleation process pre- 
dicted within this case, calculating different nucleation 
parameters. In this section, we present results for both 
metastable regions, above and below the critical chemical 
potential /icrit = 305.03 MeV, which correspond respec- 
tively to the nucleation of quark droplets in a hadronic 
environment and the formation of hadronic bubbles in a 
partonic medium. 




FIG. 1: Critical radius as a function of the quark chemical po- 
tential in the two metastable regions (between the spinodals). 
The arrows indicate the critical chemical potential. 

In Fig. [1] the critical radius, namely the radius of 
the critical bubble, is displayed. Any bubble created in 
the system via external influences or thermal fluctuations 
will either grow or shrink unless its radius equals the 
critical value, which settles the threshold between sup- 
pressed and favored bubbles. For radii bigger than the 
critical radius, the minimization of energy implies that 
the droplet will grow and eventually complete the phase 
conversion. The critical radius goes to zero at the spin- 
odal curves, where the metastable false vacuum becomes 
unstable (the barrier disappears) and the phase conver- 
sion occurs explosively via the spinodal decomposition 



process^. On the critical line the vacua become degener- 
ate so that both of them are stable yielding no nucleation, 
or equivalently, a divergent critical radius. In the case of 
the chiral transition in the LSM with quarks, we obtain 
that the critical radius is Rc < 10 fm over about ~ 90% 




FIG. 2: Ratio between the correlation length ^ and the critical 
radius Rc as a function of the quark chemical potential in the 
two metastable regions. 

The correlation length ^ which provides a measure of 
the size of the bubble wall is plotted in Fig. [21 in units 
of the critical radius. As discussed in the previous sec- 
tion, the thin-wall approximation relies on the assump- 
tion that the ratio between the surface correlation length 
^ and the radius of the critical bubble is small. Fig. [2] 
shows clearly that this is a reasonable condition in the 
vicinity of the critical line, away from the spinodals. 

The surface tension is a key parameter in quantify- 
ing the nucleation process in a given medium since it 
measures the amount of energy per unit of area that is 
spent in the construction of a surface between the phases. 
In Fig. [3] we show our result for the surface tension, 
E, for the chiral transition in the cold and dense LSM 
with quarks. The surface tension assumes values between 
^ 4 and ~ 13 MeV/fm^, being, throughout the whole 
metastable region, of the order of magnitude that renders 
the formation of quark matter viable during a supernova 
explosion, according to Ref. Q. The biggest values oc- 
cur near criticality, since this domain is characterized by 
large barriers and a small free energy difference between 
the true and false vacua. It should be noted that those 
values are well bellow previous estimates of this quantity, 
as will be discussed in the conclusion. 

We have also estimated the nucleation rate for homo- 
geneous nucleation ^] as F ~ Tje~^''/'^f , where Ft 



^ In the thin- wall approximation, which is a poor description of the 
regions close to the spinodals, the critical radius does not vanish, 
only becomes very small. The same is true for other quantities, 
as the surface tension. 
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FIG. 3: Surface tension as a function of the quark chemical 
potential in the two metastable regions. 



is the free energy of the critical bubble configuration and 
Tf = 30 MeV is an ad hoc temperature. Applying the ex- 
pression above for the nucleation rate in this case means 
that we are neglecting the temperature dependence of the 
critical-bubble free energy in the exponent. The results 
are shown in Fig. 2] for both metastable regions. The 
nucleation rate falls abruptly as the chemical potential 
approaches its critical value. 

It is interesting to point out that the difference in the 
size of the metastable regions above and below criticality 
might play an important role in determining whether nu- 
cleation is a viable process of phase conversion for a given 
system, as compared to spinodal decomposition scenar- 
ios or even if the lifetime of the system represents enough 
time for nucleation to take place. It is clear in Fig. 21 for 
instance, that the domain in chemical potential for which 
the nucleation rate assumes sizable values is larger at the 
metastable region corresponding to bubble formation of 
hadronic matter in a partonic medium. 




FIG. 4: Nucleation rate as a function of the quark chemical 
potential in the two metastable regions. 



C. Effects from vacuum terms versus thermal 
corrections 

Having presented in the subsection above the impor- 
tant parameters for nucleation and discussed their results 
for the zero-temperature chiral transition in the LSMq, 
let us now analyze the influence of quantum logarithmic 
terms in the vacuum effective potential as well as thermal 
corrections. 

In this section we focus on the metastable region above 
the critical chemical potential, aiming for thermal nucle- 
ation of quark matter droplets in compact objects and 
especially supernovae explosions. In this vein, we also 
keep relatively low temperatures when including thermal 
corrections, up to T = 30 MeV. The systematic inclusion 
of thermal corrections provides a measure of the validity 
of the zero-temperature approximation for the thermody- 
namics and phase structure in low-temperature scenarios 
such as those found in compact objects. We find that 
results for the nucleation parameters up to T = 10 MeV 
present variations within ~ 10% of the zero-temperature 
values, displayed in the previous section. 
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FIG. 5: Density in units of the nuclear saturation density 
no = 0.16 fm"'^ as a function of the quark chemical potential 
in the cold and dense LSM (curve (c)), including vacuum 
terms (curve (6)) and thermal corrections (curve (a)). 

In Figures [BHSl the role played by temperature correc- 
tions and quantum vacuum effects is presented. More 
specifically, we compare results for the three following 
situations: (a) classical vacuum effective potential with 
thermal corrections for T = 30 MeV, (6) quantum vac- 
uum effective potential at zero temperature, (c) classi- 
cal vacuum effective potential at zero temperature, i.e. 
the results shown in the last subsection. Each case 
is associated with a different critical chemical poten- 
tial (indicated by arrows in the plots; the normalization 
/^crit,o = 305.03 MeV corresponds to the critical chemical 
potential of the case (c), of the last subsection) and also 
a different spinodal chemical potential (denoted by the 
vertical, dashed lines). 
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FIG. 6: Critical radius in tiie nucleation region of quark 
droplets as a function of the quark chemical potential in 
the cold and dense LSM (curve (c)), including vacuum terms 
(curve (6)) and thermal corrections (curve (a)). 



Concerning the modification of the metastable regions, 
we show that the quantum vacuum corrections increase 
significantly the domain in chemical potential in which 
metastability persists; the metastable region in the pres- 
ence of quantum vacuum terms is 40% bigger than the 
one associated with the classical vacuum effective poten- 
tial. This feature indicates that such corrections might 
generate large differences in the dynamic evolution of the 
system, even though the absolute value of the critical 
chemical potential itself shifts only ^ 2%. 

During the dynamical process of phase conversion, the 
density is increased (e.g. via gravitational pressure in 
the supernova explosion scenario) and therefore the free 
energy required for including a new quark in the system 
also augments. The mapping between density and quark 
chemical potential within the cases we consider is plot- 
ted in Figure [5] The discontinuities signal the first-order 
phase transition that restores chiral symmetry at high 
energies. As the evolving system reaches a high enough 
density, it penetrates the metastable region above the 
critical chemical potential, where higher densities are as- 
sociated with higher supercompression. 

It is clear that the values of density predicted below 
the transition are irrealistically low. This is a conse- 
quence of the fact that the LSMq fails to describe nu- 
clear matter. Being an effetive theory constructed with 
the purpose of describing the chiral symmetry features 
of strongly interacting matter, the LSMq does not pre- 
dict the low-energy nuclear matter properties nor the gas- 
liquid transition that forms it. However, since our ulti- 
mate interest is investigating the formation of quark mat- 
ter in astrophysical processes involving ultra-compact 
objects, the framework most suited is exactly one that 
describes the high-energy transitions of strong interac- 
tions, namely chiral restoration and/or deconfinement. 
We choose therefore the LSMq and analyze the nucle- 
ation of chirally-symmetric droplets in a strongly inter- 



acting medium. 

In order to interpret the results for the nucleation pa- 
rameters, one should keep in mind the variation of the 
metastable regions and its implication on the value of 
chemical potential that needs to be reached above the 
respective criticality, i.e. the amount of supercompres- 
sion. 

Figure[6]shows the variation of the radius of the critical 
bubble, while Figure [7] displays the correlation length (or 
equivalently the approximate width of the bubble wall) 
for the three cases. A direct consequence of the shift- 
ing of the metastable domain is the fact that, for fixed 
chemical potential, the critical radius decreases signifi- 
cantly both when one includes quantum vacuum correc- 
tions (curve (&)) or in the presence of thermal effects 
(curve (a)). However, the amount of supercompression 
needed for reaching a certain viable value of critical ra- 
dius clearly increases when one includes quantum correc- 
tions and decreases in the case with thermal corrections. 
Therefore, there is a nontrivial competition between those 
effects: while quark matter nucleation is disfavored by 
the inclusion of quantum vacuum corrections, it is fa- 
vored by thermal corrections to the nucleation parame- 
ters. 
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FIG. 7: Ratio between the correlation length ^ and the critical 
radius Rc as a function of the quark chemical potential in 
the cold and dense LSM (curve (c)), including vacuum terms 
(curve (6)) and thermal corrections (curve (a)). 



The competition between vacuum and thermal correc- 
tions is more explicitly shown by our findings for the 
surface tension, displayed in Figure [H As compared to 
the zero-temperature, classical computation (curve (c)), 
the establishment of a interface between quark matter 
droplets and the hadronic medium costs more with the 
incorporation of quantum corrections in the effective po- 
tential. On the other hand, finite temperature terms tend 
to push the surface tension down, reaching a minimum 
of only ~ 3 MeV/fm^, facilitating the surface formation. 
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D. Consequences for quark- matter-induced 
supernova explosion scenario 

To contribute to the investigation of the possibihty of 
nucleating quark matter droplets during the early post- 
bounce stage of core collapse supernovae, we follow Ref. 
and define the nucleation time as being the time it 
takes for the nucleation of a single critical bubble inside 
a volume of Ikm'^, which is typical of the core of a pro- 
toneutron star, i.e. 



'^nucl 



1 



Ikm-^ 



(25) 



FIG. 8: Surface tension in the nucleation region of quark 
droplets as a function of the quark chemical potential in the 
cold and dense LSM (curve (c)), including vacuum terms 
(curve (fe)) and thermal corrections (curve (a)). 



Therefore, if the environment in the core of a compact 
object is cold enough and if the effective theory including 
quantum corrections in the vacuum is the most suited for 
describing the chiral transition, then the surface tension 
can be as high as ~ 17.5 MeV/fm^ near criticality ren- 
dering quark matter nucleation much slower and even an 
improbable phenomenon. 
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FIG. 9: Nucleation rate in the nucleation region of quark 
droplets as a function of the quark chemical potential in the 
cold and dense LSM (curve (c)), including vacuum terms 
(curve (6)) and thermal corrections (curve (a)). 

Analogous results for the nucleation rate F can be seen 
in Figure[31 Once again, the amount of supercompression 
associated with a given chemical potential should be kept 
in mind. The nucleation rate falls abruptly for chemical 
potentials approaching the respective critical values. In 
the vicinity of the spinodals, the nucleation rate is dom- 
inated by the pre-exponential factor and reaches sizable 
values, around ~ 0.1 MeV/fm'^. 



FIG. 10: Nucleation time defined for a typical volume of 1 
km'^ in the cold and dense LSM (curve (c)), including vacuum 
terms (curve (6)) and thermal corrections (curve (a)). 

Fig. [lU]shows this quantity as a function of the normal- 
ized chemical potential for the three scenarios considered 
in this analysis. The relevant time scale to compare is the 
time interval the system takes from the critical chemical 
potential to the spinodal during the supernova event, if 
it ever reaches such high densities in practice. Implicit 
in the definition above is the approximation of constant 
density and temperature over the core, which is fine for 
an estimate, since density profiled are quite flat in this 
region of the star. The typical time scale for the early 
post-bounce phase is of the order of a fraction of a second, 
so that the time within the metastable region is smaller, 
in the ballpark of milliseconds. 

Keeping these simplifications in mind, together with 
the caveat that the LSMq does not describe nuclear mat- 
ter besides the chiral properties of strong interactions, 
we show in Fig. 9 that vacuum corrections tend to in- 
crease the density depth required for efficient nucleation 
whereas thermal corrections push in the opposite direc- 
tion, consistently with the results discussed in the pre- 
vious subsection. Our results for the nucleation time, 
together with those for the surface tension, tend to favor 
the best scenario for nucleation of quark matter in the 
supernova explosion scenario considered in Ref. @, es- 
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pecially when thermal corrections with physical temper- 
atures are included. To a great measure, this happens 
because the nucleation time, and the whole process of 
phase conversion, depends very strongly on the surface 
tension, since it enters cubed in the Boltzmann exponen- 
tial of the rate F. 



IV. CONCLUSIONS AND OUTLOOK 

In this paper we have computed the effective potential 
for the linear sigma model with quarks in the MS scheme 
at zero or low temperature and finite quark chemical po- 
tential, including vacuum and medium fluctuations. We 
have discussed the issues related to the vacuum contri- 
butions and the parameter fixing in the presence of these 
corrections. Having the full effective potential, we stud- 
ied homogeneous nucleation in a framework that allowed 
for analytic calculations, given by a potential fit by a 
quartic polynomial and the thin-wall approximation. All 
the relevant quantities were computed as functions of 
the chemical potential (or the baryonic density) , the key 
function being the surface tension. 

The value of the surface tension for the QCD phase 
transitions could actually play important roles in dif- 
ferent physical phenomena. As discussed above, if this 
value is small enough, quark-matter formation could oc- 
cur during core-collapse supernovae explosions, providing 
an alternative dynamics and even an observable signal of 
the QCD phase transition within compact objects. More- 
over, a reasonably small surface tension between partonic 
and hadronic phases could also contribute to allow for dif- 
ferent compact star structures including mixed phases, as 
argued in Ref. j51|] . 

Estimates for the surface tension between a quark 
phase and hadron matter were considered previously in 
different contexts. In a study of the minimal interface 
between a color-flavor locked phase and nuclear matter 
in a first order transition, the authors of Ref. fs^ use di- 
mensional analysis and obtain S ^ 300 MeV/fm^ assum- 
ing that the transition occurs within a fermi in thickness. 
Taking into account the effects from charge screening and 
structured mixed phases, the authors of Ref. [53] provide 
estimates in the range of 50 — 150 MeV/fm^ but do not 
exclude smaller or larger values. 

In this paper we computed the surface tension as a 



function of quark chemical potential (or as a function of 
baryon density) within the linear sigma model, isolating 
the role played by quantum vacuum terms and thermal 
corrections. In particular, we show that the model pre- 
dicts a surface tension of S ~ 5-15 MeV/fm^, rendering 
nucleation of quark matter possible during the early post- 
bounce stage of core collapse supernovae. Including tem- 
perature effects and vacuum logarithmic corrections, we 
find a clear competition between these features in char- 
acterizing the dynamics of the chiral phase conversion, so 
that if the temperature is low enough the consistent in- 
clusion of vacuum corrections could help preventing the 
nucleation of quark matter during the early post-bounce 
stage of core collapse supernovae. 

As discussed in Subsection III El effects from interac- 
tions between the sigma field and the quarks that come 
about at two-loop order could in principle contribute 
to this competition as a third sizable modification and 
should be investigated as well. 

The linear sigma model, however, does not contain es- 
sential ingredients to describe nuclear matter, e.g. it does 
not reproduce features such as the saturation density and 
the binding energy. Therefore, the results obtained in 
this paper should be considered with caution when ap- 
plied to compact stars or the early universe. It is an 
effective theory for a first-order chiral phase transition in 
cold and dense strongly interacting matter, and allows 
for a clean calculation of the physical quantities that 
are relevant for homogeneous nucleation in the process 
of phase conversion. In the spirit of an effective model 
description, our results should be viewed as estimates 
that indicate that the surface tension is reasonably low 
and falls with baryon density, as one increases the su- 
percompression. First-principle calculations in QCD in 
this domain are probably out of reach in the near future. 
Therefore, estimates within other effective models would 
be very welcome. 
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